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Objective Cholinergic nicotinic receptor (CHRN) gene family has been known to mediate the highly additive effects of nicotine in the 
body, and implicated nicotine dependence (ND) and related phenotypes. Previous studies have found that CHRNA6-CHRNB3 cluster 
polymorphisms were significantly associated with the risk of ND and various tobacco behaviors. The aim of study was to evaluate the 
genetic association of CHRNB3 and CHRNA6 polymorphisms with the risk of ND based on the Fagerstrom Test for Nicotine Depen- 
dence (FTND) score and five subscales of nicotine dependence syndrome scale (NDSS) in Korean population. 

Methods Six SNPs in CHRNA6-CHRNB3 cluster were analyzed in 576 Korean subjects. Association analysis using logistic models and 
regression analysis with NDSS were performed. 

Results There was no association in the case-control analysis, whereas all six SNPs were significantly associated with drive factor among 
NDSS in subgroup based on the FTND score. CHRNB3 rs49S4 and CHRNA6 rsl6891604 showed significant associations with NDSSF1 
(drive) in dominant models among moderate to severe ND among smokers after correction (p c ""=0.02 and 0.001, respectively), whereas 
other four SNPs showed significant associations among mild ND after correction (p'°"=0.03-0.02 in dominant model). 

Conclusion This study showed that the genetic influence of CHRNB3-CHRNA6 cluster polymorphisms are found in a ND endophe- 
notype (drive) using NDSS subscales, rather than the risk of ND in Korean population. Our findings might be the first report for the as- 
sociation of CHRNB3-CHRNA6 cluster with ND-related phenotypes in Korean and might offer an approach to elucidating the molecu- 
lar mechanisms of ND and ND-related phenotypes. Psychiatry Investig 2014;11(3):307-312 
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INTRODUCTION 

Cigarette smoking is a major health problem and a leading 
environmental risk factor for several diseases, such as pulmo- 
nary disease and cancer worldwide. Despite the increasing 
awareness of health-related consequences, a large number of 
smokers fail to stop smoking due to their nicotine dependence 
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(ND). 1 The genetic influences on nicotine dependence were 
established through twin studies, 2 4 and genetic factors with 
heritability account for between 40-70% of the phenotypic 
variance. 3 " 6 

Cholinergic nicotinic receptor (CHRN) gene family code 
pentameric molecular assemblies of nicotinic acetylcholine 
receptor (nAChR) subunits mediate the highly additive effects 
of nicotine in the body. 7 According to the functional impor- 
tance of nicotinic subunit genes, early genetic association stud- 
ies of nicotine dependence have been focused on CHRNA4 
and CHRNB2, 7,S and the significant associations of the CHR- 
NA5-CHRNA3-CHRNB4 cluster, CHRNB3-CHRNA6 cluster 
and level of nicotine dependence have been recently report- 
ed. 910 Besides these subunit genes, CHRNB3 and CHRNA5 
have emerged as strong candidates for nicotine dependence 
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in two recent reports of a genome-wide association and high- 
throughput candidate gene survey, 910 and variants in or near 
CHRND-CHRNG, CHRNA7, and CHRNA10 show modest 
association with nicotine dependence risk in African-Ameri- 
cans. 11 It has been also suggested that genetic variants in the 
CHRN gene family are associated with tobacco -related behav- 
iors. CHRNA5 rsl6969968 and rs588765 affect smoking quan- 
tity, 12 and CHRNA4 variants have been associated with tobac- 
co-related behaviors. 81314 CHRNB3 rs4950 and rsl3280604 
are significantly associated with the subjective response factors 
to initial tobacco use. 15 Although studies on smoking and ND 
based on linkage scans and the candidate gene approach have 
yielded somewhat inconsistent results, 16 genetic influences of 
the CHRN gene family have been replicated in separate com- 
munity samples and other ethnic populations. 11,15 

The assessment of ND has relied largely on the use of the 
Fagerstrom Tolerance Questionnaire (FTQ) or the Fagerstrom 
Test for Nicotine Dependence (FTND). 17,18 The Fagerstrom 
scales were originally developed as unidimensional measures 
of physical tolerance. There is no assessment of other impor- 
tant aspects of ND, such as craving, subjective compulsion to 
smoke, nicotine withdrawal, behavioral saliency, or behavior- 
al automaticity, which are now regarded as core constructs. 
The nicotine dependence syndrome scale (NDSS) was devel- 
oped using a multidimensional concept of ND, and includes 
scores on five subscales: Drive, Priority, Tolerance, Continu- 
ity, and Stereotypy. 19 

Previous studies have found that the CHRNA6-CHRNB3 
cluster polymorphisms were significantly associated with the 
risk of ND and various tobacco behaviors, 9,10,15,20 and these as- 
sociations were replicated in a large meta-analysis in samples 
of European ancestry from the ENGAGE consortium. 21 To the 
best of our knowledge, there is no genetic association study of 
the CHRNA6- CHRNB3 cluster with the risk of ND and the 
multidimensional nicotine dependence measure in the Asian 

Table 1 . Clinical profiles of study subjects 



population, especially in the Korean population. The aim of 
the current study was to evaluate the genetic association of 
CHRNB3 and CHRNA6 polymorphisms with the risk of ND 
based on the FTND score and five subscales of NDSS in the 
Korean population. 

METHODS 
Subjects 

Nicotine-dependent patients in several community mental 
health centers and hospital nicotine dependence clinics were 
recruited. All (n=371) were diagnosed using the Diagnostic 
and Statistical Manual of Mental Disorders, Fourth Edition 
(DSM-IV) criteria. The population controls (n=205) consist- 
ed of randomly selected students from Eulji University and 
Ulsan University of Korea. After a detailed explanation about 
the process of this study, informed consent was obtained. 
The Institutional Review Board of the hospital had approved 
the study. Clinical data from the K-NDSS tests (Korean ver- 
sion of NDSS) were obtained. The clinical profiles of the study 
subjects were summarized in Table 1. 

SNP selection and genotyping 

Candidate SNPs were selected using the database in the 
Asian population from the International HapMap Project da- 
tabase (http://hapmap.ncbi.nlm.nih.gov). SNP selection was 
based on the following criteria; In the 80-kb region including 
the CHRNB3 and CHRNA6 genes, SNPs were selected taking 
into consideration their minor allele frequencies (MAFs; >5%) 
and linkage disequilibrium (LD) (only one SNP if there were 
in absolute; r 2 >0.95). 

For genotyping of the polymorphic sites, the amplifying 
primers and probes were designed for TaqMan (Supplementa- 
ry Table 1, online). The PCRs were run on a TaqMan Univer- 
sal Master Mix without UNG (Applied Biosystems, Foster 



Smoker 



Description 


A11ND 


Severe+Moderate D 
(FTND>4) 


MildND 
(FTND<4) 


NS 


Number of subjects 


371 


183 


188 


205 


Age (min-max) 


44.16(20-72) 


46.44 (20-70) 


41.94 (20-72) 


20.94(18-37) 


Gender (M/F) 


362/9 


181/2 


181/7 


102/103 


NDSSF1 (Drive) 


13.53±3.65* 


14.73±3.19 


12.26+3.68 




NDSSF2 (Priority) 


5.88±2.92 


6.36±3.07 


5.35+2.65 




NDSSF3 (Tolerance) 


6.41 ±2.06 


6.1512.16 


6.69+1.92 




NDSSF4 (Continuity) 


5.51±2.25 


6.0312.23 


4.95+2.15 




NDSSF5 (Stereotypy) 


5.97±2.02 


6.48+1.95 


5.42+1.96 





*meanistandard deviation. ND: nicotine dependent subject, NS: non-smoker of population control, FTND: Fagerstrom Test for Nicotine 
Dependence, NDSS: Nicotine Dependence Syndrome Scale 
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City, CA, USA) but with a PCR primer concentration of 900 
nM and a TaqMan MGB-probe concentration of 200 nM. 
The reactions were carried out in a 384-well format in a total 
reaction volume of 5 uL with 20 ng of the genomic DNA. The 
plates then were placed in a thermal cycler (PE 9700, Applied 
Biosystems) and heated to 50°C for 2 minutes and 95°C for 10 
minutes; this step was followed by 40 cycles of 95°C for 15 sec- 
onds and 60°C for 1 minute. Ten percent of the samples were 
randomly selected for repeated genotyping, and the results 
were 100% concordant. The TaqMan assay plates were then 
transferred to a Prism 7900HT instrument (Applied Biosys- 
tems), which we used to measure the fluorescence intensity in 
each well of the plate. The fluorescence data files from each 
plate were analyzed with automated software (SDS 2.1). 

Statistical analysis 

Haplotypes were inferred using the algorithm developed 
by the Broad Institute, Haploview. 22 We examined Lewontins 
D' (|D'|) and the linkage disequilibrium coefficient r 2 be- 
tween all pairs of biallelic loci. 23 The allelic distributions of 
polymorphisms and haplotypes among patients with nicotine 
dependence and normal subjects were evaluated by logistic 
regression models calculating odd ratios (OR), 95% confi- 
dential intervals (CI) and corresponding P-values using age 
and sex as covariates. According to the FTND scores, the gen- 
otype distributions of patients with severe and mild nicotine 
dependence were also compared to that of non-smokers. A re- 
gression model was used for association analyses with five 
NDSS subscale factors, adjusting for age and sex as covariates. 
To correct for multiple testing, the effective number of inde- 
pendent markers was calculated using the software SNPSpD 
(http://genepi.qimr.edu.au/general/daleN/SNPSpD), which 
is based on the spectral decomposition (SpD) of matrices of 
pair-wise LD values between SNPs. 24 The resulting number 
of independent marker loci of the CHRNB3-CHRNA6 region 
(3.499) was applied to correct for multiple testing. We ana- 
lyzed all the subjects and the dichotomized subgroups; mild 
ND (FTND less than 4) and moderate-to-severe (FTND given 



4 or more). 

RESULTS 

In this study, we examined the association of polymorphisms 
in the CHRNB3 and CHRNA6 gene region with the risk of 
nicotine dependence in Korean subjects (n=576). Asian gen- 
otype information on 72 SNPs in about 80 kb regions of the 
CHRNB3 and CHRNA6 gene were extracted from the Inter- 
national HapMap database, and six SNPs were finally selected, 
and then genotyped in 371 smokers and 205 non-smokers. 
The minor allele frequency and Hardy- Weinberg equilibrium 
are summarized in the Supplementary Table 2 (online). Link- 
age disequilibria among SNPs were measured by calculating 
Lewontins D and r 2 values (Supplementary Figure 1, online), 
and rsl3280604, rsl689161, rs9298628 and rsl0109040 were 
tightly linked with one another. The haplotypes were con- 
structed using six SNPs and only one haplotype (htl) among 
three major haplotypes was used for further analyses because 
ht2 and ht3 were almost equivalent to rsl6891604 and rs4954, 
respectively. 

Initially, when the association analysis of six SNPs and a 
haplotype were examined using logistic regression models be- 
tween smokers and non-smokers, no association was observed 
(Table 2). The genetic effects of polymorphisms on the five 
subscale factors of NDSS were analyzed using regression mod- 
els. CHRNA3 rs4954 and CHRNA6 rsl6891604 showed sig- 
nificant associations with NDSSF1 (drive) in dominant mod- 
els among moderate-to-severe ND among smokers after 
correction (p m "=0.02 and 0.001, respectively) (Table 3), whereas 
the other four SNPs were significantly associated with NDSSF1 
among mild ND after correction (/>' D "=0.04-0.03 in the co- 
dominant model and p" 5n =0.03-0.02 in the dominant model) 
(Table 3). Overall associations showed that the drive subscale 
of ND was decreased in subjects with minor alleles, but 
CHRNB3 rs4954 showed an increase in the drive subscale. In 
the additional association analyses for the rest of the NDSSFs, 
no association was observed (Supplementary Table 3, online). 



Table 2. Association analysis of CHRNB3 and CHRNA6 polymorphisms with the nicotine dependence in Korean subjects 



Gene 



Loci 



Allele 



NDS (N=371) NS (N=205) 



OR (95%CI) 



Statistical power (%) 



CHRNB3 


rsl3280604 


A>G 


0.173 


0.210 


0.94 (0.34-2.63) 


0.91 


58.66 




rsl6891561 


OT 


0.171 


0.205 


0.98 (0.34-2.83) 


0.97 


57.95 




rs4954 


A>G 


0.115 


0.141 


1.08 (0.30-3.92) 


0.90 


46.70 


CHRNA6 


rs9298628 


OT 


0.190 


0.227 


1.01 (0.35-2.97) 


0.98 


60.88 




rs 16891604 


C>A 


0.218 


0.176 


1.13 (0.45-2.82) 


0.80 


53.40 




rs 10109040 


C>G 


0.182 


0.212 


0.83 (0.27-2.53) 


0.74 


58.93 




htl 


Ahtl>htl 


0.588 


0.590 


0.94 (0.42-2.12) 


0.88 


72.95 



Odd ratio and corresponding p-value were calculated by using logistic models whilst controlling by age and sex. NDS: nicotine dependent 
smoker, NS: non-smoker of population control, CHRN: cholinergic nicotinic receptor, OR: odd ratio, CI: confidence interval 
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Table 3. Regression analyses of CHRNB3 and CHRNA6 polymorphisms with NDSSF1 (drive) among smokers 



Group Gene 


Loci 


C/C 


C/R 


R/R 


Pa 


Pa co " 


Pb 




Pc 


Pc" 1 " 


A11ND CHRNB3 


rsl3280604 


238 (13.6113.55) 


96 (13.33+3.98) 


12 (13.67+3.06) 


0.35 


NS 


0.35 


NS 


0.70 


NS 




rsl6891561 


239(13.6213.55) 


95 (13.31+3.99) 


12 (13.67+3.06) 


0.33 


NS 


0.32 


NS 


0.70 


NS 




rs4954 


270(13.60+3.51) 


69 (13.14+4.24) 


7(14.86+2.54) 


0.81 


NS 


0.70 


NS 


0.73 


NS 


CHRNA6 


rs9298628 


226(13.7113.53) 


107 (13.0813.96) 


13 (14.08+2.96) 


0.14 


0.50 


0.09 


0.30 


0.98 


NS 




rs 16891604 


205 (13.71+3.92) 


127 (13.2913.26) 


13(13.31+2.93) 


0.15 


NS 


0.10 


0.34 


0.97 


NS 




rs 10109040 


232(13.71+3.53) 


100 (13.0813.98) 


14 (13.79+3.04) 


0.14 


0.49 


0.11 


0.39 


0.68 


NS 




htl 


67(13.21+3.27) 


155 (13.34+3.64) 


124 (13.94+3.85) 0.02* 


0.07 


0.10 


0.36 


0.03* 


0.10 


Severe & CHRNB3 


rsl3280604 


122 (14.42+3.28) 


49 (15.47+3.01) 


8 (15.13+2.30) 


0.16 


NS 


0.11 


0.38 


0.82 


NS 


moderate ND 


rsl6891561 


123 (14.43+3.27) 


48 (15.4613.04) 


8 (15.13+2.30) 


0.18 


NS 


0.12 


0.43 


0.82 


NS 




rs4954 


143 (14.40+3.20) 


30 (16.27+2.82) 


6(15.17+2.64) 


0.02* 


0.07 


0.006* 


0.02* 


0.76 


NS 


CHRNA6 


rs9298628 


114(14.6113.25) 


56 (14.9313.22) 


9(15.11+2.15) 


0.59 


NS 


0.60 


NS 


0.79 


NS 




rs 16891604 


102(15.29+3.31) 


72 (13.94+2.95) 


5 (14.80+1.30) 


0.004* 


0.02* 


0.001* 


0.005* 


0.99 


NS 




rs 10109040 


118(14.6213.29) 


51 (15.04+3.09) 


10 (14.60+2.59) 


0.67 


NS 


0.51 


NS 


0.80 


NS 




htl 


40 (14.48+2.58) 


78 (14.4613.27) 


61 (15.26+3.42) 


0.07 


0.25 


0.38 


NS 


0.04* 


0.15 


Mild ND CHRNB3 


rsl3280604 


116 (12.75+3.64) 


47 (11.11+3.65) 


4 (10.75+2.22) 


0.008* 


0.03* 


0.007* 


0.03* 


0.37 


NS 




rsl6891561 


116(12.7513.64) 


47(11.11+3.65) 


4 (10.75+2.22) 


0.008* 


0.03* 


0.007* 


0.03* 


0.37 


NS 




rs4954 


127 (12.6913.63) 


39 (10.74+3.54) 


1 (13.00) 


0.009* 


0.03* 


0.006* 


0.02* 


0.84 


NS 


CHRNA6 


rs9298628 


112(12.79+3.57) 


51 (11.06+3.72) 


4(11.7513.50) 


0.01* 


0.03* 


0.005* 


0.02* 


0.72 


NS 




rsl6891604 


103(12.14+3.86) 


55 (12.44+3.46) 


8 (12.3813.34) 


0.64 


NS 


0.57 


NS 


0.99 


NS 




rs 10109040 


114(12.77+3.54) 


49(11.04+3.79) 


4(11.7513.50) 


0.01* 


0.04* 


0.007* 


0.02* 


0.72 


NS 




htl 


27(11.33+3.33) 


77 (12.21+3.67) 


63(12.67+3.82) 


0.13 


0.46 


0.16 


0.58 


0.26 


0.91 



C/C, C/R and R/R represent major homozygotes, heterozygotes and minor homozygotes, respectively. Pa, Pb and Pc mean the p-values of co- 
dominant, dominant and recessive model, respectively. *p<0.05. ND: nicotine dependence, NS: non-smoker of population control, CHRN: cho- 
linergic nicotinic receptor, NDSS: Nicotine Dependence Syndrome Scale 



DISCUSSION 

Drive, which is an essential factor in addiction and one of 
the NDSS factors, is defined as a presence of craving or with- 
drawal in the abstinence period. Using the regression analy- 
ses of subscales of NDSS among Korean smokers revealed that 
CHRNB3 and CHRNA6 polymorphisms were significantly 
associated with the drive subscale factor (Table 3). Interesting- 
ly, the genetic associations of four SNPs in CHRNB3 and CHR- 
NA6 {rsl3280604, rsl6891651 in CHRNB3; and rs9298628, 
rs 10109040 in CHRNA6) were found only in mild ND, where- 
as CHRNB3 rs4954 and CHRNA6 rsl6891604 showed the as- 
sociations in moderate-to-severe ND. The genetic effects in the 
dominant model of CHRNB3 rs4954 and CHRNA6 rsl6891604 
were opposite from each other, that is, CHRNB3 rs4954 in- 
creased the craving of smoking in moderate-to-severe ND 
subjects but CHRNA6 rs 1689 1604 reduced it. Our results 
showed that none of the CHRNB3 and CHRNA6 polymor- 
phisms affected the diagnosis of ND; however, CHRNB3 
rs4954 and CHRNA6 rs 1689 1604 might play an important 
role in the craving and withdrawal of smoking in the absti- 
nence period of severe ND subjects, respectively. On the other 



hand, the rest 4 SNPs of CHRNB3 and CHRNA6, which sh- 
owed tight linkage with one another, showed the genetic ef- 
fect for the withdrawal of smoking. These results implicated 
that the genetic function of CHRNB3-CHRNA6 gene cluster 
polymorphisms might differ in accordance with the phase of 
ND. Although CHRNB3 rs4954 also showed the reduced ef- 
fect of the drive subscale in mild ND, and this effect of CHRNB3 
rs4954 might be due to its low number of subjects with minor 
homozygote. The limitations of the current study were the 
low number of subjects and the differences in the male/female 
ratio and average age between ND smokers and non-smokers. 
Although these limitations could affect our results by positive 
or negative false errors, our results might offer a plausible clue 
of the genetic factors related with ND-related endopheno- 
types. In fact, the statistical power in the case-control analysis 
was relatively low (46.70-72.95%) (Table 2), and the additional 
analysis of the risk of ND in male subjects showed that there 
was no big difference in the distribution and association re- 
sults (Supplementary Table 4, online). 

Several studies have shown that polymorphisms in the CH- 
RNB3-CHRNA6 cluster are associated with the risk of ND and 
related phenotypes, such as smoking behaviors. CHRNB3 
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rs9298629 (absolute LD with rs9298628 in the current study) 
showed nominal association with the risk of ND in female Is- 
raelis. 20 Associations of SNPs in the 5' region of CHRNB3 and 
3' UTR of CHRNA6 with the risk of ND were also found in 
US-based samples and Australian-based European- Ancestry 
samples. 9 Moreover, CHRNB3 rs4950 and rsl3280604 were 
significantly associated with ND-related phenotypes, such as 
the subjective responses to nicotine. 15 However, the association 
of previous studies showed somewhat different results to our 
results, especially the association results for the risk of ND. 
The analysis of the MAF and LD structure of 6 polymorphisms 
among three typical ethnic groups, including Asian, Caucasian 
and African data from 1000 genome databases (http://br- 
owser.1000genomes.org/index.html) revealed that the MAF 
of CHRNA6 rs 1689 1604 was quite different from that of 
Caucasians and Africans (Supplementary Table 5, online), and 
that the LD structures of tested ethnic groups were quite dif- 
ferent with one another (Supplementary Figure 2, online). Th- 
erefore, our results could be valuable because polymorphisms 
tested might be specific to ethnic group. 

In a study of NDSS using a criterion-group design contrast- 
ing chippers and regular smokers, the drive subscale was the 
most robust correlate of group membership, with an odds ra- 
tio of 6.45, accounting for 54% of the variance. The area under 
the ROC curve for Drive was 0.92, again indicating a very high 
degree of discrimination. 25 Therefore, the drive subscale may 
be the most important factor among the subscales of NDSS. 
Exploring the appropriate phenotype or endophenotype in a 
genetic study is essential; however, it is difficult to know in 
advance which phenotype might be associated with specific 
genes or SNPs. The use of multidimensional scales such as 
NDSS potentially aids in exploring the construct and nature of 
ND. In a recent Finnish cohort study, NDSS were used as an 
ND endophenotype measurement, and it was reported that 
the tolerance factor of the NDSS showed a suggestive associ- 
ation to CHRNB4 rsl 1636753, rs 11634351, and rsl948. 26 

The nicotine dependence phenotype/endophenotype could 
be various and differ by gender. As mentioned in the limita- 
tions, our study did not have an even male/female ratio. In fu- 
ture research, age- and sex-matched genetic analysis is needed. 

An a priori assumption of our study was the genetic influ- 
ences of CHRNB3 and CHRNA6 for the risk of ND as well as 
ND-related endophenotype. Our results implicated that CH- 
RNB3 and CHRNA6 polymorphisms might be genetic pre- 
dispositions for the ND-related endophenotypes rather than 
the risk of ND in the Korean population. To the best of our 
knowledge, our results may be the first report investigating the 
association of the CHRNB3-CHRNA6 cluster with ND-relat- 
ed endophenotypes, including subscale factors of NDSS in 
the Korean population. The current study may be meaningful, 



although there were several limitations in our study (e.g., the 
small sample size and smaller portion of females among ND 
subjects than that of controls). 

Conclusion 

In conclusion, this study showed that the genetic influence 
of CHRNB3-CHRNA6 cluster polymorphisms are found in 
the ND endophenotype (drive) using NDSS subscales, rather 
than the risk of ND in the Korean population. These effects 
of multiple loci were different between mild and moderate- 
to-severe groups of ND. Our findings could offer an approach 
to elucidating the molecular mechanisms of ND and ND-re- 
lated phenotypes in the Korean population. 
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